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Abstract 

In simulations of (magnetized-) fluid dynamics in physics and astrophysics, the visu- 
alization techniques are so frequently applied to analyse data that they have become 
a fundamental part of the research. Data produced is often a multi-dimensional set 
with several physical quantities, that are usually complex to manage and anal- 
yse. JETGET is a visualization and analysis tool we developed for accessing data 
stored in Hierarchical Data Format (HDF) and ASCII files. Although JETGET has 
been optimized to handle data output from jet simulations using the Zeus code from 
NCSA, it is also capable of analysing other data output from simulations using other 
codes. JETGET can select variables from the data files, render both two- and three- 
dimensional graphics and analyse and plot important physical quantities. Graphics 
can be saved in encapsulated Postscript, JPEG, VRML or saved into an MPEG for 
later visualization and/or presentations. An example of use of JETGET in analysing 
a 3-dimensional simulation of jets emanating from accretion disks surrounding a 
protostar is shown. The strength of JETGET in extracting the physics underly- 
ing such phenomena is demonstrated as well as its capabilities in visualizing the 
3-dimensional features of the simulated magneto-hydrodynamic jets. The JETGET 
tool is written in Interactive Data Language (IDL) and uses a graphical user inter- 
face to manipulate the data. The tool was developed on a LINUX platform and can 
be run on any platform that supports IDL. JETGET can be downloaded (including 
more information about its utilities) from http://www.capca.ucalgary.ca/software. 
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Complete paper (including high resolution figures) can be found at 
http://www.capca.ucalgary.ca/software/jetget/ 

Fig. 1. The Main JETGET module. This is the main JETGET window. Prom here, 
all of JETGETs modules are accessible. 

1 Introduction 

Today the scientific community is confronted with the problem of understand- 
ing and analysing large masses of numeric scientific data thanks to the avail- 
ability of super-computing systems. One solution to this problem is scien- 
tific visualization: converting the numeric data into pictures more readily un- 
derstandable by scientists. However, while visualization helps understand the 
overall behavior of the simulations, it is still a daunting task to extract the 
underlying physics. Ideally one would like to use a tool where both visualiza- 
tion and analysis are built in. Particularly, in many cases in computational 
(magnetized- )fiuid dynamics in physics and astrophysics, extensive data pro- 
cessing is required to obtain meaningful information. 

JETGET is a data visuahzation and analysis software specifically designed to 
deal with data from such simulations. JETGET is IDL (Interactive Data Lan- 
guage) based. It provides a set of modules we built to aid the scientist when 
analysing and visualizing 2-dimensional (2-D and 2.5-D) and 3-dimensional 
(3-D) data. It can handle large datasets allowing both their graphical repre- 
sentation and analysis. JETGET allows users to interact with time histories, 
profiles, contour and surface plots on the same window, and quickly perform 
mathematical manipulations or combinations of data. 

Here we describe the basic functionalities of JETGET. Note that a help file 
for each module is provided which contains more detailed information. This 
paper is organized as follows. We start in section 2 by a description of the 
Setup modules where specific information about the data files is specified. In 
sections 3-10 the different modules of JETGET arc described. In particular, we 
show how the modules can be used to extract various and crucial information 
from the simulations. A discussion of the different output methods is given in 
section 11 before concluding in section 12. 



2 Setup and Data Input 

The Main module in JETGET is shown in Fig. 1. From this window, the differ- 
ent modules of JETGET are accessible. In this paper, we use a low resolution 
(60 X 55 X 55 zones) version of the simulation performed in [1] to illustrate the 
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Complete paper (including high resolution figures) can be found at 
http://www.capca.ucalgary.ca/software/jetget/ 

Fig. 2. The Setup module. The figure shows the main Setup module. From Setup, 
two other modules can be accessed, to define the structure of the data files, and to 
define new variables. 

many functionalities of these modules. The simulations in question ^ describe 
a jet emanating from the surface of an accretion disk surrounding a protostar. 
The plot on the main widget shows the central star and the twisted magnetic 
field configuration at later stages of jet evolution from this simulation. Carte- 
sian coordinates (xi, X2, X3) were used for this simulation; the corresponding 
indices are (i,j,k). The disk is taken to be in the X2-Xs (or j-k) plane, and the 
disk axis corresponds to the xi-axis (or i direction; see Fig. 1 in [1]). Variables 
are in dimensionless units, as defined by the problem to be simulated (e.g. 

[I])- 

The Setup module (Fig. 2) is the configuration tool for JETGET. Here the 
general information concerning the data files for JETGET is set. First the 
user must specify the path to the directory where JETGET and the data-files 
are stored. Then, the dimension (two or three) of the data-files must be given, 
since this affects the structure of the data files (see section 2.2). JETGET has 
originally been designed to handle 3-D HDF files. However, few modules such 
as Scan/Slices, can also handle 2-D files. Thus it is necessary to specify (in 
the Setup module) the dimension of the dataset. 

JETGETs default nomenclature is set according to data output from ZEUS 
simulations with the default names "hdfaaOOOOOO.xxx" (with x's being digits 
defining the file number). JETGET expects the file names to be of the form 
"AAAxxxBBB" , where "xxx" indicates the file number. The number of digits 
in the filename is given to JETGET in the Setup module. The "AAA" and 
"BBB" (called prefix and postfix in the Setup module) are characters, also 
given to JETGET in the Setup module. These fields can be left blank, so that 
the numbers can come first or last in the filename (dots are also allowed). The 
time-step (dt) between successive files is also defined in the Setup module. It 
is necessary to give the file ID number for the first file in the dataset. This file 
will be read to determine the scaling and resolution of the simulations. 

When pressing the "Save" button in any of the other JETGET modules, the 
content of the render window will be saved to a Postscript (PS) or a JPEG 
file. The choice, PS or JPEG is made in the Setup module. The image will 
be saved with colors. If Postscript is chosen, the resolution of the postscript 
file is also passed from the Setup module. The resolution is given as Dot Per 
Inch (DPI). For printing there are a few important points worth mentioning. 

^ The corresponding output/HDF files analysed in this paper can be downloaded 
from jhttp:/ / www.capca.ucalgary.ca/software^ 
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Complete paper (including high resolution figures) can be found at 
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Fig. 3. The Primitive Variables module. In this module it is possible to define the 

structure of the data files, if different from the default structure. Shown is the default 
structure for 3 dimensional files. Also, it is possible to choose between HDF and 
ASCII data files. Here HDF files is selected. 

Complete paper (including high resolution figures) can be found at 
http ://www. capca.ucalgary . ca/ software/j etget/ 

Fig. 4. The Calculate Variables module. Here the user can define new variables not 
initially present in the data-files. 

The printer has a maximum resolution, hence choosing a larger resolution than 
what the printer can handle will not be beneficial. Furthermore, the size of the 
postscript file depends strongly on the resolution, so limiting the resolution is 
a good idea. Especially output with color will produce large files. The PS and 
JPEG output methods are described further in section 11. 

The Panel Setup configures the Panel which is an alternative display feature 
(see section 11.1). The number of rows and columns to use in the Panel is set 
in Setup. 

2.1 Coordinate system 

Most 3-D jet simulations are run using Cartesian coordinates, even though 
cylindrical coordinates might seem more intuitive. As discussed in [1], there 
are many challenging numerical problems when using cylindrical coordinates. 
JETGET therefore requires Cartesian coordinates for 3-D data files. However, 
it can automatically calculate the corresponding poloidal (e.g. Bp) and toroidal 
(e.g. Bff) components of the magnetic field and the velocity field. In Fig. 11, 
the toroidal component of the velocity {V^) is shown in two slices through the 
dataset. In Fig 12 B^/Bp is plotted along a magnetic field line. 

2.2 Primitive variables: default and user variables 

The structure of the data files (Fig. 3) is given in the Primitive variables 
module accessible from the Setup module. JETGET requires all variables from 
one time dump to be put in the same file. The default file format is Hierarchical 
Data Format (HDF), although ASCII ^ files can also be analysed. 

^ When using ASCII files the user must specify: (i) The size of the grid, for instance 
60x60x60 zones. This should be an array containing the two or three sizes, (ii) The 
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JETGET can handle both 2-D and 3-D data files. It can handle both strictly 
2-D files and 2.5-D files. The difference is that 2.5-D files have 3 components 
of the vector variables (e.g. the magnetic field), whereas 2-D files only have 
2 components. 2.5-D files have been used in jet simulations, assuming axial 
symmetry. 

For 2-D [3-D] datafiles the default structure is: Ist-component of the velocity 
(Vi), 2nd-component of the velocity (V2), [3rd-component of the velocity (V3)], 
Ist-component of the magnetic field (-Bi), 2nd-component of the magnetic field 
(S2), [3rd-component of the magnetic field (-B3)], density (d), specific energy 
1 (ei; e.g. gas pressure), specific energy 2 (62; e.g. turbulent pressure). 2.5-D 
datafiles have a default structure similar to 3-D files. 

When using the default structure, some additional variables are defined. In 
the 3-D case, these variables are the logarithm of the density, the logarithm of 
the two specific energy components (ei and 62), the kinetic energy, magnetic 
energy, the poloidal component of the velocity (Vp) and magnetic field (Bp), 
the toroidal component of the velocity (V^) and magnetic field {B^) and the 
inclination angle of the magnetic field with respect to the surface of the accre- 
tion disk. In the 2-D and 2.5-D case, these extra variables are the logarithm 
of the density and the two energy components. 

If the structure of the data files is changed, the correct structure must be given 
in the Primitive Variables module. The name of each variable must be given 
in the right order. If a variable is a component of a vector, the vector name 
must be given as well. 

2.3 Calculate variables 

The Calculate variables module (Fig. 4) is accessible from the Setup module. 
It allows the user to define new variables based on the primitive variables. 
This could for instance be the square of the magnetic field {B^) or the square 
of the velocity field (V^) as seen in Fig. 4. It is also possible to calculate a 
new set of variables making up a vector quantity. These new vectors will then 
be available for visualization using for instance the Vectors or the Fieldlines 
modules. 

If the Save button is clicked in the Setup module after calculating new vari- 
ables, the definitions for these variables will be stored in the ".jetget" config- 
uration file, see section 2.4. 

scaling in each direction. This is a one dimensional array with the size of the grid. 
So two or three of these arrays are needed, (iii) The actual data-arrays, in the order 
given in JETGET's Primitive Variables module. 
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Complete paper (including high resolution figures) can be found at 
http://www.capca.ucalgary.ca/software/jetget/ 

Fig. 5. The Fieldlines module. The image shows 50 magnetic field lines anchored 

to a rotating Keplerian disk. The field lines twist in the central region. In this case 

the footpoints of all the field lines are on a line. This, as well as the number of field 

lines drawn, can be changed in the Fieldlines Setup module, accessible from the 

Fieldlines module. 

Complete paper (including high resolution figures) can be found at 

http://www.capca.ucalgary.ca/software/jetget/ 

Fig. 6. The Fieldlines Setup module. This widget configures the initial set up for 
the visualization of the field lines in the Fieldlines module. Shown (as crosses) are 
the footpoints or starting points in the disk of the 50 magnetic field lines shown in 
Fig. 5. 

2.4 Storing information 

In each of the different JETGET modules, there is a "Save" button. In the 
Setup module pressing this button will save the setup in the ".jetgct" config- 
uration file in the directory specified by the environment variable HOME (i.e. 
the home directory on a UNIX computer) and is automatically read next time 
JETGET is started. 

Information given to the Fieldlines module can be stored in a configuration 
file named ". jet get -fieldlines" , also located in the directory specified by the 
environment variable HOME. 



3 Fieldlines 

The Fieldlines module (Fig. 5) will draw field lines of any vector field in the 
dataset. These will be visualized in three dimensions. Therefore, this module 
will not work for a 2-D dataset. Fig. 5 shows the visualization of 50 selected 
magnetic field lines anchored to the accretion disk. The footpoints are deter- 
mined by the Fieldlines Setup module (section 3.1. The field lines are found by 
using first order forward stepping and linear interpolation in three dimensions. 

Since the field lines are found using a first order technique, there is some 
uncertainty on the calculated field lines. The user can check how accurately 
the plotted field lines have been computed by pressing the Check button. This 
will compute the field lines, first by stepping outwards from the disk, and then 
by taking the endpoint of the computed field lines and step back towards the 
disk. The two sets of field lines will then be drawn side by side, for qualitative 
inspection. Also, the difference in the footpoints for the field fines found by 
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stepping outwards from the disk and end points for the field fines found by 
stepping backwards towards the disk wiU be given, as weU as the maximum 
difference and the mean difference. 

The agreement of the field lines when stepping from the disk out, and stepping 
back towards the disk depends on the complexity of the field lines being visu- 
alized. Tests calculating magnetic field lines from the simulations performed 
in [1] have shown that there is a good agreement provided the field lines do 
not extend too far out into the coarse grid. Solutions should be checked for 
each plot using the "Check" button. If the accuracy is not high enough, a 
smaller step length in the Fieldlines Setup module (Fig. 6) should be used. 

The rendering of field lines can be used to show that these magnetic field 
lines are wound up and coUimated in simulations of jets. Similarly, velocity 
streamlines ^ can be drawn. For an example of the use of magnetic field lines 
and streamlines, see for instance Fig. 14 and 15 in [1]. 

3.1 Fieldlines Setup 

In the Fieldlines Setup module (Fig. 6), accessible from the Fieldlines module, 
the visualization properties of the field lines can be changed. The number of 
footpoints (the number of field lines to be drawn) are set in this module, as 
well as the layout of the footpoints. Fig. 6 shows 50 footpoints located on a 
line in the rotating accretion disk. 

In the Fieldlines Setup module it is also possible to specify the inner and outer 
radii for the footpoints. This allows the study of only field lines originating 
from a particular part of the disk and its influence on the outflow. As seen 
in Fig. 5 the outer fleld lines are almost straight. In some cases these may be 
considered undesirable, and by choosing a smaller radius in which to draw the 
field fines, these will be avoided. 

The viewing angle of the field lines can be selected in the Fieldlines Setup 
module, for an optimum view of the field lines. The view can be defined with 
respect to any of the three major axes. We note that the field lines can not be 
rotated once drawn. To see the field fines from another angle, they need to be 
regenerated. Fig. 6 shows how the starting points will be placed in the view 
with the current configuration, together with the disk axis perpendicular to 
the footpoints. This image is drawn whenever the "Use" button is pressed. 

The configuration information given in this module can be stored in a file in 
the home directory called ".jetget_fieldlines" by clicking the "Save" button. 

^ With streamline we mean the field line of the velocity field. 
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Complete paper (including high resolution figures) can be found at 
http://www.capca.ucalgary.ca/software/jetget/ 

Fig. 7. The Fluxes module. The image shows the Mass flux through slice 30 in the 
data set in the jet (xi) direction. 

Complete paper (including high resolution figures) can be found at 
http://www.capca.ucalgary.ca/software/jetget/ 

Fig. 8. The Fourier module. The image shows the kr vs m contours at six diff'erent 
times. Notice how the m = mode is dominant in the earlier stage. Later, the 
171 = 1 mode appears. 

This file is automatically read every time the Fieldlines module is started. 



3. 2 Limitations 



The footpoints of the field lines to be drawn must be anchored in the xi = 
plane (the accretion disk in the simulations analysed here; see Fig. 6). 



4 Fluxes 



The Fluxes module (Fig. 7) will plot the flux vs time through a slice in the data 
set perpendicular to the disk axis. It can plot the mass flux, the momentum 
flux, the energy flux and the magnetic flux. The slice is chosen as a slice 
number in the data file. 

JETGETs Fluxes module provides a direct calculation of the fluxes from jet 
simulations, and these can then be directly compared to observed fluxes. The 
authors in Ref. [2] and [3] have calculated and plotted the mass flux, the 
momentum flux and the kinetic energy flux in their 2.5-D simulations of non 
relativistic outflows from Keplerian accretion disks. Such plots can now easily 
be generated by using JETGETs Fluxes module. Fig. 7 shows the mass flux 
through slice 30 in the simulations analysed here. It can be seen from this 
flgure, that the flow reaches this slice roughly at time 50. At time 150 a 
decrease in the mass flux is noticeable due to the slowing down and the sideway 
motion of the jet. 
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Fig. 9. The Fourier module showing mode amplitudes. Negative m values are shown 
as a dotted line, whereas positive m values are shown as a solid line. The figure shows 
amplitudes for m = ±5. 

5 Fourier 



The Fourier module (Fig. 8) calculates the dominant MHD modes as the 
simulation evolves in time. 

Given a field F(r) (here pressure) on a Cartesian grid x=(a;2,a;3) where r = 
\jx\^ Xg, the Fourier transform gives 

1 

!m{r) = F{r) X exp{-im(p)d(f), 

ZTT Jo 

and it is possible to find the m-modes, given by: 

fm{r)w{r)rdr 

Jo 

where w{r) is the filter. It consists on the decomposition (Fourier transform) 
of the two dimensional pressure distribution of the selected slices. 

To run the Fourier module, one must first define a filter for the deconvolution. 
The filter depends on the ro parameter, which is the jet's average radius and on 

the dvf) parameter, which is the deconvolution window. While there are two 
pre-defined filters, a Sine-square and a Gaussian, other filters can be easily 
added. 



5.1 kr — m contours 

When the first two steps of the Fourier module have been performed, one 
can visualize the kr — m contours, kr is the radial wave number. The authors 
in Ref. [1] used this module to extract the dominant MHD modes in their 
jet simulations. This showed the dominant MHD modes (in particular the 
m = 1 kink mode) in the jet evolution and the importance of mode coupling 
in stabilizing the jet. 
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Complete paper (including high resolution figures) can be found at 
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Fig. 10. The Isosurface module. The image shows the isosurface of Vi, taken at a 
level of Fi = 0.1 at time = 80. 

5.2 Show mode amplitudes 

The given mode amplitude can be plotted for one or more files. Negative m 
value is shown as dotted line, whereas positive m value is shown as a solid 

line, as seen in Fig. 9. The mode is given as an input parameter to the Fourier 
module. The mode shown is the relative amplitude a{m)/a{m = 0). 



6 Isosurface 

The Isosurface module (Fig. 10) will draw isosurfaces of a selected quantity. 
The viewing angle can be set freely, however, once a viewing angle is chosen 
and the object drawn, this can not be rotated without regenerating the object. 
The evolution of a given surface can be followed over time, either by making 
consecutive plots shown one at a time, or several plots shown side by side in 
a panel display. 

In Fig. 10 an isosurface of Vi (the velocity component along the disk axis) 
is shown together with the Isosurface widget. The isosurface has been drawn 
with filled colors and shading. Another option is to draw using a wired mesh 
generating a more transparent isosurface. 

6. 1 Limitations 

The view is always a square, so aspect ratios are not maintained for different 
viewing angles. 



7 Mosaic 

The Mosaic (Fig. 11) module will draw filled contour plots of slices through a 
2-D or 3-D dataset. This module will draw three variables at the same time. 
In a 3-D dataset, the Mosaic module will show a slice along the jet, and 
one perpendicular to the disk axis. In a 2-D dataset, the entire slice will be 
shown, still for three different variables. All plots will be scaled relative to 
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Complete paper (including high resolution figures) can be found at 
http://www.capca.ucalgary.ca/software/jetget/ 

Fig. 11. The Mosaic module. The image shows slices taken both along the disk axis 
(left panel) and perpendicular to the disk axis (right panel). The three variables 
plotted are F^, Bi and the logarithm of the density. The cut along the disk axis 
is taken at slice 27, the cut perpendicular to the disk axis is at slice 29; in the 
middle of the simulated box. The jet assumes an elliptical shape at this stage in 
the simulation, indicating the importance of the m = ±2 MHD modes, as already 
indicated in Fig. 8. 

Complete paper (including high resolution figures) can be found at 
http://www.capca.ucalgary.ca/software/jetget/ 

Fig. 12. The Quantities module showing B^p/Bp vs distance along field line. The 
footpoints are (x2,X3) = (3,2). For S < 40, it can be seen that B^f, dominates Bp, 
whereas for S > 40 Bp is dominating B^. The narrow peak at about 5 = 40 is 
typical for the Alfven surface. 

the physical scale of the grid. This means that the size of the window will be 
adjusted accordingly. 

In this module it is not possible to select contour plots, only filled plots. But 
as in the Scan/Slice module, there is a wide range of color palettes to choose 
from. The contour levels will be determined from the maximum and minimum 
levels for each sequence of shces in each panel. 

Fig. 11 shows a slice along the jet and a slice across the jet for each of the 
three variables V^, Bi and Log{d). The Xi-x^ slice along the jet correspond to 
j — 29, the X2-Xs sUce across the jet is taken eik i — 27, in the middle of the 
simulated box. 



8 Quantities 

The Quantities module (Fig. 12) will plot any selected quantity along a field 
line of any vector quantity in the data set. To find the field line, a footpoint 
in the disk must be provided. This footpoint is given in physical coordinates. 
The step length and the number of steps to be used when finding the field line 
to plot the quantity along must also be given. The algorithm for finding the 
field fine is the same as for the Fieldlines module (see section 3) . 

In [4] the possibility that energy and angular momentum are removed mag- 
netically from accretion disks are examined. They studied the evolution of 
different quantities along a magnetic field line. This can now easily be ex- 
tracted from numerical simulations using the Quantities module. The image 
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Complete paper (including high resolution figures) can be found at 
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Fig. 13. The Scan/Slice module. The figure shows plots of the density at different 

times using a panel as output. The slices is taken parallel to the disk, at slice 17. 
One can clearly see that the jet twists and become deformed, and in the later images 
it is about to regain its circular shape. 

in Fig. 12 show a plot of the B^/Bp along a magnetic field line starting at 
(2^2, 2:3) = (3.0,2.0); very close to the innermost radius of the disk. The loca- 
tion of the Alfven surface (where the jet speed equals the Alfven speed) at 
5" = 40 is evident in this plot. Within the = 40 radius the jet morphology 
remains cylindrical and axi-symmetric with the m = and m = ±2 modes 
clearly dominant. Beyond the Alfven surface, the m = ±1 (kink) mode takes 
over reducing the strength of the toroidal magnetic field component with re- 
spect to the poloidal component. 

8. 1 Limitations 

The Quantities module uses the same engine to find field lines as the Fieldlines 
module, thus the same limitations apply to the Quantities module as to the 
Fieldlines module. In particular, field lines evolving too far into a coarse grid 
might not be accurate. 



9 Scan/Slice 

The Scan/Slice module (Fig. 13) displays contour plots or filled plots of a 
selected variable. The variable can be scanned in either time (the same slice 
in different files) or space (different slices in the same file) . The contour levels 
will be determined from the maximum and minimum levels of all the selected 
slices. This means that if the variable changes rapidly throughout the dataset, 
some plots will only show a few contours. This way it is easier to compare 
the plots. If a constant number of contours were drawn for each plot, small 
features would be over-exposed. 

If filled plots are chosen, a built-in IDL routine named XPalette can be used 
to select and manipulate the colors. The user can choose between a wide range 
of built-in color tables in IDL. A color bar will also be drawn. 

The Scan/Slice module does not need a uniform grid. In the outer region 
of a simulation, a coarse grid (a grid with increasing zone sizes) will often 
be introduced to increase the physical size of the grid without increasing the 
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Fig. 14. The Vectors module. The image shows the velocity vectors portraying the 
rotation of the jet, taken in a slice parallel to the disk at slice 18. Only high velocity 
vectors have been drawn. 

number of grid points too much. Scan/Slice will take care of the proper scaling 
of the plot. 

The Scan/Slice module handles both two and three dimensional data files. 
This way of visuahzing data files has been used many times, see for instance 
[2] where 2.5-D MHD simulations of the evolution of non relativistic outfiows 
from Keplcrian accretion disks orbiting low mass protostars or black holes 
accreting at sub-Eddington rates arc treated. Or [1] where three dimensional 
MHD simulations of the outflow from magnetized accretion disks are treated. 



10 Vectors 

The Vectors module (Fig. 14) is used to visualize the vector quantities in the 
dataset. The visualization will show arrows representing the two components 
of the vectors in a shce/plane [{x2,Xs) in Fig. 14]. The length of the vectors 
can be adjusted by a scale factor, since the vectors often will be too small to 
see. Therefore, the length of the vectors should only be viewed as relative to 
other vectors, not as a measure of the absolute length of the vectors. 

It is possible to limit the plotted vectors to a certain range of lengths. This 
is convenient when for instance only plotting long vectors, thus getting rid 
of small irrelevant vectors. Note that this selection is made before the length 
scale factor is applied to the vectors. Furthermore, it is possible to plot only 
a subset of the vectors in the chosen range. If a vector is drawn in each point, 
the drawn vectors might overlap and make the plot difficult to see. Several 
plots can be made in sequence, either one at a time or several side by side in 
a Panel display. The vectors can be scanned in either time (the same slice in 
different files) or space (different slices in the same file). 

While the Fieldlines module is efficient for visuahzing the vector field in 3-D, 
the Vectors module will show the magnitude of the fields, but only in a 2-D 
slice. The image in Fig. 14 shows velocity vectors taken in a slice perpendicular 
to the disk axis. The image clearly shows the clockwise rotation of the jet 
material. Examples of the use of two dimensional vector plots can be seen in 

[!]• 
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Complete paper (including high resolution figures) can be found at 
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Fig. 15. The Filename input pop up window. This window will pop up when a file 
is to be saved. In this case an eps file is to be saved, therefore the choice of vector 
or bitmap postscript is available. 

11 Output 

There are several different output methods from JETGET. Common for all 
modules are the render window, where the requested output will be drawn 
when clicking the "Run" button. (In the Fourier module there is no run button, 
since there are two different visualization methods, hence two buttons labeled 
"Show k-m contours" and "Show mode amphtudes" exists.). If more than one 
object is to be drawn, the second object will overwrite the first one. 

When the "Save" button in one of JETGETs modules is pressed, a window 
(Fig. 15) will pop up asking for the file name to save to. Note that file name 
extensions (e.g. .jpg) are not automatically added by JETGET. 



11.1 Panel 



The Panel is an alternative display method for some of the modules and allows 
for several plots to be visualized simultaneously. The panel will draw a set of 
different plots side by side, rather than one at a time. Information for the size 
of the Panel is passed from the Setup module, see section 2 (e.g.. Fig. 13). 



11.2 Postscript (PS) 

If PS is chosen in the Setup module, the content of the render window will be 
saved to an encapsulated Postscript file. The quality of the PS file (the dpi; the 
number of dots per inch) can be given in the Setup module, as discussed in sec- 
tion 2. The Postscript files can be written either as vectors or as bitmaps. The 
difference being that vectors can effectively be resized later, whereas bitmaps 
cannot. However, it is important to note that vector Postscript files in some 
cases become huge (50 MB). The choice (vector or bitmap) is made in the 
Filename input pop-up window (Fig. 15) as where the filename is given. 
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11.3 Joint Photographic Experts Group (JPEG) 

If JPEG is selected as the output file type in the Setup module, the content 
of the render window will be saved to a JPEG file when the save button is 
pressed. 

11.4 Moving Picture Experts Group (MPEG) 

MPEG output is another possibility in the modules where several objects can 
be drawn consecutively. The views are then stored in an MPEG file. The 
generated MPEG file will be made with 30 frames per second, and each frame 
is stored 30 times in the file, so that each image will be visible for one second 
when the MPEG movie is played. 

11.5 Virtual Reality Modeling Language (VRML) 

VRML output is possible from the Fieldlines module and from the Isosurface 
module. VRML is a file format allowing 3-D objects to be viewed on the 
world wide web (with the appropriate plug-in). When viewing a VRML file, 
the object can be rotated, and one can zoom in and out on the object. Only 
one object can be put into one VRML file. 



12 Conclusion 

JETGET can be used to analyse and visualize most data output from simu- 
lations of (magneto-)hydrodynamic fiuids, specifically simulations of jets. Be- 
sides the reconstruction of the three dimensional shapes and features of the 
jets, JETGET allows us to extract vital information concerning their ener- 
getics, dynamics and stability. This leads to a deeper understanding of jet 
physics. JETGET is free for download and can be found at 
|http: / /www.capca.ucalgary.ca/software, . 
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A Requirements 
A.l IDL 

JETGET uses newly implemented calls in IDL, thus IDL 5.6 or newer is 

required. If the MPEG feature in several of the modules should be used, a 
special MPEG license is necessary. This must be obtained separately from 
RSI. IDL runs on several UNIX systems, Linux systems, Windows and MAC 
systems. JETGET was developed and has been tested on IDL 5.6 for Linux. 

A. 2 Memory us age /requirement 

IDL itself takes up about 10 MB of memory. The memory needed for each 
module depends strongly on the size of the data files. Using the panel option 
with filled contours in the Scan/Slices module is probably the module requiring 
the most memory. With a 60 x 60 x 60 data file, this can amount to 30 MB. 
For a 200 x 200 x 200 data file, more than 100 MB is needed. 

A. 3 Color depth 

The recommended color depth is 16 bit. Other depths will work, but there are 
problems with backing store. That is, that if an other window is placed on top 
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of the render window in JETGET, the contents of the render window should 
be redrawn when the other window is moved. 
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